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Qh' Abstract 
(D 

' Gauge theories induced by scalars in the fundamental representation of the U {Nc)ga.ugc x U (A''j)giobai 

group are investigated in the large Nc and Nf limit. A master field is defined from bilinears of the 
scalar field following an Eguchi-Kawai type reduction of spacetime. The density function for the 
^ I master field satisfies an integral equation that can be solved exactly in two dimensions {D = 2) 

' and in a convergent series of approximations at D > 2. While at D = 2 the system is in the same 

^ . phase at all e = Nc/Nj, it undergoes a phase transition at a critical value, ec{D), for D > 2. 

a^ ■ 
o 

^ ! 1 Introduction 



Oh! 



X 



Many quantum field theories, with the notable exception of nonabelian gauge theories, can be 
solved exactly when the number of fields, A^, tends to infinity. The large N expansion is not overly 
useful for solving nonabelian gauge theories because the leading order of the expansion is equivalent 
. to the sum of all planar diagrams. Unfortunately, it is not much simpler to sum planar diagrams 

than to solve the exact theory. 

To obtain exact results further simplifications, like those proposed by Kazakov and Migdal Q, 
are required. Using the lattice action of a gauged adjoint representation scalar field, in which the 
5^ I pure gauge term was omitted, these authors derived an equation for the master field, composed of 

the eigenvalues of the scalar field matrix. Once the master field is known other physical parameters 
of the system can also be calculated. Following the philosophy that nonabelian gauge theories are 
unique one hopes that correct self-interactions of the gauge field are generated by scalar dynamics. 

There are two reasons why gauge theories induced by adjoint representation scalars |^ 
may not be completely equivalent to standard nonabelian theories. The first reason is well-known: 
Because the action is quadratic in the compact gauge field, U, there is a hidden local Z^- symmetry, 
which implies that the vacuum expectation value of Wilson loops is zero. Q 

The second reason which makes one believe that induced theories and standard theories are 
different, goes as follows. Suppose that along with Araf 'eva Q we gauge only one of the indices of 
the adjoint representation scalar field. Then the theory has a U{Nc)ga,uge x f^(-^c)giobai symmetry 
and the unwanted symmetry is not present. We can view such a theory as one having global 
and gauge symmetries with equal numbers of fiavors and colors. The gauge theory induced by 
such a scalar field is not asymptotically free in four dimensions and as such is not likely to be 
confining either. In fact, the induced theory will only be asymptotically free if the number of 
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complex scalars, Nj, satisfies the inequafity Nj > llA'c This is the primary reason why we 

previously considered a scalar field theory with U{Nc) gauge x f^(-^/)giobai symmetry. Q 

The prevalent argument against theories with fundamental representation scalars, rather then 
adjoint representation ones has been the lack of existence of a master field, precluding chances 
for finding exact solutions. While in gauge theories induced by adjoint representation scalars the 
master field is composed of the eigenvalues of the scalar field matrix, for fundamental representation 
scalars no such quantity seems to exist. 

The purpose of this paper is to show that the general bias against gauge theories induced 
by fundamental representation scalars is unfounded. It is possible to introduce an adjoint gauge 
tensor formed from bilinears of the fundamental field that can play a role identical to that of the 
primary scalar fields in previously investigated induced theories. A master field is formed from the 
eigenvalues of this positive definite self-adjoint matrix. The density function for the components of 
the master field satisfies a singular integral equation that can, in particular cases, be solved exactly 
or in a convergent series of approximations. 



2 Lattice model for gauged fundamental scalars 

Following the philosophy of induced gauge theories the Wilson plaquette term is omitted from a lat- 
tice action in the hope that correct gauge interactions are generated though scalar self- interactions. 
Then the only term of the Lagrangian in which the gauge field makes its appearance is the gauged 
kinetic term. The kinetic terms written in terms of fundamental representation scalars, ipl^{x), and 
compact gauge fields, U^{x), has the form 

D 

Lk = hY1 [4i^)Ui.ix)'ilj'ix + ae^) + 4ix + ae^)Ulix)ij\x)'\ , (1) 
11=1 

where k is the hopping parameter, a is the lattice constant, i is a flavor index, while color indices 
are suppressed. D is the number of dimensions and is a unit vector along the /ith coordinate 
axis. 

The hopping parameter can be eliminated by rescaling the scalar field and absorbing the change 
into the coefficients of the scalar potential. For that reason in what follows we will set n = 1. 

The Lagrangian also contains interaction terms. We only introduce -0^ interaction terms but 
other terms relevant or marginal in less than four dimensions could also be introduced and treated 
without additional difficulty. In principle there are two fourth order self interaction terms allowed 
by the U{Nc) x U{Nf) symmetry, differing in their symmetry structure. These terms are formed 
from pairs of bilinears of the fundamental field that are global scalars and gauge tensors or scalars, 
or alternatively formed from gauge scalars that are global tensors or scalars. They have the form 

Lj = Lt + Ls = ^4{xW {x)^]{x)^\x) + ^^\{x)^\x)^Ij]{xW{x). (2) 

Terms with scalar coupling are not generated in the hopping parameter expansion. Thus, it is 
self-consistent to keep the term with the tensor structure only. Therefore, in what follows we will 
omit the scalar term and use the notation qt = 9- 

The complete Lagrangian of the induced theory has then the form: L = Lk + Lt + Lm, where 
Lm is the mass term, 

LM = m^4{xW{x). (3) 



2 



Following Araf'eva Q we use an Eguchi-Kawai type realization Q of coordinate space in the 
internal symmetry space. Such a reduction is admissible when the number of colors and flavors 
is large. In the large A'^^ and Nc limit one can realize the spacetime translation operators on 
the symmetry groups as D^[k) and A^(p), where D^{k) is a diagonal U{Nc) matrix, with matrix 
elements e'^"'^'^ , 1 < a < Nc, and A^(p) is a diagonal U{Nf) matrix, with matrix elements e*'*^'', 
1 < i < Nf. These operators translate objects in the fundamental representation of the SU{Nc) 
[SU (Nf)] group by one lattice unit along the /ith axis. The model thus obtained called "quenched" 
because the integration over the momentum variables should be performed after physical amplitudes 
are calculated. Then the lagrangian is reduced to a single site form 

L = n?4i^' + J-V^IVV]^' + E Tr [i:^^{p)4D]/\k)U^Dy\kW + h.c] , (4) 

where, for the sake of symmetry, we assume that the gauge fields occupy sites at the middle of 
links, halfway between two vertices of the lattice, where the scalar particles reside. It can be shown 
that the perturbation expansion of this model coincides with that of the original model at large Nf 
and Nc- 

At this point it can be easily seen that the translation operators can be eliminated from the 
theory altogether making the theory effectively translation invariant. The operator A^(A;) can be 
eliminated, because in gauge invariant closed loops its contributions cancel. The operators D^{k) 
can be eliminated using a redefinition of the gauge field. Then we arrive at the lagrangian 

L = m^^tl;' + ^4^^^]^' + ^ [ipjU^i^' + h.c] , (5) 

which reproduces the correct perturbation expansion. 

The theory defined by Lagrangian (^) fails in the weak coupling limit because it does not 
reproduce the correct eigenvalue spectrum of the matrices U^. Therefore, one should introduce, 
along with Gross and Kitazawa a quenched theory, containing a gauge invariant modification of 
the measure constraining the eigenvalues of U^. The appropriately modified theory will be discussed 
briefly at the end of this paper and in more details in a future publication. Throughout this paper 
we will concentrate on the simple theory defined by Lagrangian (||). 

At this point there is still no obvious candidate for a master field. Notice, however that la- 
grangian (^) depends only on the positive self adjoint tensor formed from bilinears of the scalar 
field, ^> 

= CV'f , (6) 

where, for the sake of clarity the gauge indices a and /? are displayed. 

We now change variables by introducing the following factor into the partition function 

1 = J d^6{^ - V-Vl) = J d^dpe^'P^^-^'^^^\ (7) 

Integrating over the complex scalar fields we obtain 



At large Nf we can evaluate the the integral over p using the saddle point method. Thus, to terms 
of 0(1/ Nf) we obtain 

"dV' ~ / (i«>e^/T'-'°s*. (9) 
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We are now in the position of writing down the complete action in terms of the fields $ and 
U. Action S will be defined with a sign opposite to the usual one so that 5'=logZ, where Z is the 
Boltzmann factor. We obtain 



S = Tr 



iV^ logc^ - m2 (D - -|- $2 + ^ (D ([/^ + C/t) 
iV f — 



(10) 



The partition function is invariant under a unitary transformation of the positive hermitian 
matrix, that diagonalizes The action, S, written in terms of the eigenvalues of 

a = l,...,iVc is 



EM 



2 + Nf\og(t)a -m?(, 



+ DW{(t>), 



where 



VF(0) = log J dUe 



$(C/+C/t) 



(11) 



(12) 



is the single link integral studied by Brower, Nauenberg, Brezin and Gross (Pfl^ in the large A'^ 
limit, and D is the number of dimensions. The first term of the action is the logarithm of the 
Jacobian of the transformation from <I> to its eigenvalues, (pa- 

As it will soon become clear the action is of 0{Nj). Then the most probable eigenvalues can 
be calculated by using the saddle point method. The values of (f)a that maximize the action form 
the master field. By definition, all eigenvalues, 4>a, are positive. Zero eigenvalues are excluded by 
the term log<I> of the action. At large N these eigenvalues approximate a continuous distribution. 
As usual, we introduce the density of eigenvalues, pi^i'), by the relation 



1 



Nc 



where v is the running eigenvalue. Consequently, p{i') satisfies the sum rule 

A^. 



p{v)du 



N 



f 



(13) 



(14) 



The range of v is yet unknown. The positivity of matrix <I> requires that it starts at some fmin > 0. 

After rescaling variable by the relation v — > I'Nf, Nf can be eliminated and the action will 
have a nontrivial dependence on e only. Using the results of |^ [|l^] and neglecting trivial constant 
terms we can write the action as 



S f [1 /" 

— = y dup{u) -j^ j d\ p{\) \og{v - Xf +\ogu - w?v - gv^ 



+ DW, 



where 



W 



dvp{u) I 2{iy^ + c 



,1/2 



dX />(A)log (1/2 + c)i/2 + (A2 + c) 



a/2 



The constant c is determined from the equation 



dv'- 



2 + c)l/2 



2. 



(15) 



(16) 



(17) 
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if condition 



du^-—^ > 2. (18) 
u 

is satisfied. This condition defines the strong coupUng phase. If s < 2 then c = and the system 
is in the weak couphng phase of the single site integral. 

An integral equation for the optimal density distribution is obtained if one varies ( |T5|) with 
respect to the distribution function 

J dnp{fi) log{p -vf - D \og{\/u'^ + C+ ^ p"^ + c) + log 1/ - rn^u + 2D\/u'^ + c-gv^ = \, (19) 



where A is a variational parameter, corresponding to constraint (|1^ 

The derivative of ( p!^ ) with respect to v provides a singular integral equation that one hopes to 
solve for p{x), 

P / dMf^) + _ _ ^2 ^ 2D 2gu = 0. 20) 



The solution of integral equation ([2^) does not determine the range of eigenvalues (variable u). 
In fact, one can get a solution for all choices of < t'mm < i^max- We need to impose an additional 
constraint, p{i') > 0, to get a unique solution. The positivity of /o(i^), as the density of eigenvalues 
is fairly obvious. Formally, it is required by the uniqueness of eigenvalues (pa of matrix $ as a 
function of a. 

In principle we could obtain equations for Vmin and Vmax by varying the action with respect 
to these parameters. Being dependent on dup{u) only, these variations lead to (^), taken at the 
upper or lower limits and multiplied by the value of plv) taken at the same points. It follows 
then that the variation of the action has a double zero at a density p{x) satisfying the variational 
equation and also vanishing at the end points. Consequently, the action has a point of inflection 
as a function of the end point parameters. If one slightly varies fmm (or Vmax) from the points 
at which p{i'rnin) = pi^max) = 0, with intention to increase S, then the density function becomes 
negative near the end point, leading to an inadmissible solution. Consequently, one must demand 
the vanishing of the density function p(i/) at the end points. As we will see later this condition 
determines the end points uniquely. 

The method we use for determining the range of variation of eigenvalues can be applied to a 
wide range of problems involving matrix variables in the large N limit and master fields. If the 
action is of the form 



k i=l 



Fk{uu...,Uk), (21) 



then the spectral function must vanish at the boundaries, Vmin and Vmax- In turn, this requirement 
can be used to determine Vmin and fmax- 

We have not been able to find an exact solution of (|2^) at arbitrary D, either in the weak or in 
the strong coupling region. At D = 2 an exact solution can be found in the weak coupling phase 
(c = 0). Furthermore, at D > 2 the integral equation can be transformed into a Fredholm type 
equation providing a convergent expansion for p{i'), displaying critical behavior, around D = 2 and 
a, probably continuous, phase transition as a function of e. The solution of the Fredholm equation 
at > 2 is also possible using a convergent expansion in a power series of e. 
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3 Exact solution of the model at D = 2 



If the number of dimensions, D = 2, then the problem of finding the density function simpHfies 
considerably in the weak coupling phase of W{a). To see this let us consider ( pO[ ) in the weak 
coupling phase 



P J dfipifi) 



D 



^ 1^ + fj- 

f2 _ m^~2D 



+ ^-M'^ -2u = 0, (22) 



where M is the renormalized mass defined by M = ™^ — . By redefining the mass parameter we 
scale out the coupling constant, g, from action. In the weak coupling phase it only appears in the 
rescaled relation (|l^) 

s = ^ J dv^^ < 2. (23) 

By an appropriate choice of g it is always possible to tune s to an arbitrary value inside the weak 
coupling phase, e.g to the boundary of the strong and weak coupling phases, where s = 2. Note 
that in the strong coupling regime c approaches zero continuously when the critical point s = 2 
is approached. Therefore, if s is tuned to the critical point then our solution is valid, no matter 
which side the critical point is approached from. 



At D = 2 the derivative of variational equation ( |22D can be transformed to 



where the constants A, z and variable x have been introduced by the relations 



A{1 + ZX), V^ax = \J A{1 + Z), Umin= \JA{1- Z),. (25) 

and where we define R{x) = p{v) In other words, the range of variable xis— l<x<l. 

Determining end points Vmax and Vmin is equivalent to determining A and z. As we will see 
later these parameters are fixed by the requirement of i?(l) = -R(— 1) = 0. 

A solution of integral equation ( p^ ) that vanishes at the end points, x = ±1, can be obtained 
by inverting the Cauchy principal value integral. After obvious manipulations one obtains the 
following expression for the density function 

^(-) = -4^^r , / +M^ + 2JA{l + zy)\ 

^JA z Vl -^2 /•! 2dC 



2 7r2 7_i (^r+i^+yr+^)v/r^ 

+ ^^^1^ r (26) 

Att'^ ^^Jl + X zJ^iyi+^Z + yJl + X Z) VI + i Z^l^ 

The density, -R(x), is obviously positive definite and as such, physically admissible. We could 
express R{x) in a closed form as a combination of complete elliptic functions. As it is not especially 
illuminating that form will not be presented in this paper. 

At this point A and z are yet undetermined. It is comparatively easy to see that consistency 
conditions allow us to determine these constants. Before discussing these conditions it is worth 
mentioning that solutions of integral equation (^), not vanishing at x = ±1, could be found at 
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every choice of the end points. These solutions would either not maximize S or would not be 
positive definite. 

To determine the end points I'min and I'max (i-e. parameters A and z) first note that form 
( p^ seems to be independent of M^. This is so because when the Cauchy integral is inverted the 
term on the right hand side of (|26| ) gives a vanishing contribution. This term can however be 
recovered if we substitute expression (p6D into (24) and use the Poincare-Bertrand formula for the 



exchange of the order of Cauchy integrals. Then all terms of the right hand side of (24) except 
are recovered. A constant term appears that must be identified with M^: 

1 /•! dC / 1 



2JA {l + z 0] = M'. (27) 



This equation can be used to determine A in terms of z and M. 

The second condition is obtained from the normalization of the spectral function, (14). After 
simple manipulations we obtain the following relation: 

1 /-i e / 1 



2JA {l + zO]=e. (28) 



2^7-1 v/T^ \VMT+YO 

Having determined A as a function of z and M from ( p7| ) we can use ( pSj ) to determine z as a 
function of the physical parameters M and e. 

Numerical calculations show that one can find unique positive solutions for A and z at all 
possible choices of real and positive e. As a contrast, we will see later at -D > 2 we can find a 
similar type of solution only in a restricted range of e. Z) = 2 is a lower critical dimension of the 
theory. 

We now examine the behavior of the solution in the limits of small and large e. At small values 
of € = Nc/Nf z ~ e. In other words, the range of eigenvalues relative to the average value shrinks 
to zero. Then the matrix $ tends to a constant multiple of the unit matrix. We can also see by 
inspection that in the limit z ^ the distribution of eigenvalues follows Wigner's semicircle law. 

2 

At large e 2 ~ 1 — 6 ^ , the lower bound of the eigenvalues of matrix ^min — (-^ — 2;) ~ 
ee~^ , approaches zero, while the upper limit, i'max = A {1 + z) 2e, tends to infinity, so that 
the eigenvalues span the whole positive real axis. The distribution of eigenvalues, R{x), differs 
substantially from the semicircle law, favoring eigenvalues close to the lower bound. 



4 Induced theories at D > 2 

At > 2 we have not been able to find an exact solution for the density function. It turns out 
however that even at D > 2 the singular integral equation can be transformed into a Fredholm 
equation, solutions of which can be generated in convergent expansions. Having these expansions 
is almost as useful as having an exact solution. 

We will start again from the derivative of (^) with respect to v. Using variable x and parameters 
A and z as at -D = 2, in ( p^ we obtain 

1 , R(x) (D-2) z R(x) dx 

dx- — ' 



1 y - X 2 J-l + Z X {y/l + Z X + ^J1 + z y) 

1 



^/_j^^ + M'^ + 2V^(l + .x), (29) 
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This equation is identical to (24) other then the appearance of the first term on the right hand 
side. 

After inverting the Cauchy integral on the left hand side of (^) we obtain 



Rix) 



f{x) + {D-2) j\yK{x,y) R{y) 
\fA z VI - a;2 /•! 2di 



2 7r2 



-1 {^l + i z+ ^l + x z)^l^ 

1 de 



4 ^2 VI Vl + 2; ^ J-i Wl + iz + .Jl + x z)s/l + i Z^/l^ 
{D - 2) z2 Vl -x2 /-i i?(y) dy 



4 7r2 7-1 (VI + y z + VI + 2; 2) VI + y ^ 
-1 (vi + g 2 + V1 + 2; ^)(Vi + C + Vi + y -2)71^ 



(30) 



This is a linear, inhomogeneous integral equation with an inhomogeneity that is identical to the 
density function at D = 2, with kernel 



K{x,y) 



4 7r2 {^l + y z + y/l + x z) yjl + y z 
1 dg 



-1 (vi + e ^ + vi + :e z)(vi + e ^ + vi + y •2)vi^' 



(31) 



and with expansion parameter D — 2. Kernel K{x, y) can easily be symmetrized. It is bounded, 
i.e. it satisfies condition 

rl rl 

(32) 



J dx J dy\K{x,y)\'^ < 00. 



It follows that (|3^) is a Fredholm equation. 

The solution of the integral equation is given by the following expression 

R{x) = fix) + r dy Res{x,y) f{y), 



(33) 



where f{x) is the inhomogeneity of the integral equation, and Res{x,y) is the resolvent of the 
equation, 

{D-2)K{x,y) N{x,y) 



Res{x, y) 



(34) 



I- [D -2) K{x,y) Det 

The Fredholm numerator, N{x,y), and the Fredholm determinant, Det, are entire functions of the 
expansion parameter, D — 2. The first terms of their expansions are given by 



iV(x,y) = iD-2)K{x,y) 



Det 



1-{D- 2)TvK + 



{D - 2)2 
2! 

P-2) 
2! 



[K{x,y) TTK-K\x,y))+... 



(35) 
(36) 



where is understood in the operator sense. 
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The negative definite kernel insures that at D > 2 the Fredholm determinant has no zeros. At 
zeros of the Fredholm determinant only the homogeneous equation has a solution. Thus, the only 
reason why solutions of (|30| ) would not be acceptable is the lack of positivity at all x. 

Note that at z = the kernel vanishes. Then, in view of the positivity of the inhomogeneity 
at all —1 < X < 1 the spectral function is physically acceptable. Since small z and small e are 
synonymous, we find that at all -D > 2 one can find a critical value of e, such that for e < €c{D) a 
physically acceptable solution exists. As we saw earlier ec(2) = oo. 

By examining the solution at D = 2 one can see that z ^ 1 corresponds to e ^ oo. Therefore, 
one would expect that at D > 2 limitations on the range of admissible values of z may appear near 
z = 1 (large e), excluding large values of e. Therefore we performed a systematic expansion of the 
Fredholm determinant and solution R{x) near z = 1. 

Due to singularities at x = —1 the terms of Det (and of N{x,y)) diverge logarithmically near 
z = 1. The Fredholm determinant can be written as 



Det = exp{Trlog(l - {D - 2) K)} = exp |-(Z) - 2) TrK - 



-TtK^ + ...]. (37) 



It is easy to calculate the leading (and even subleading) logarithms of the traces of powers of the 
kernel. Using these we obtain the following expression for the Fredholm determinant near z = 1: 



Det = {I - zY^'^\ (38) 

where 

1 rn-9 ('n-9'i2 

(39) 



a{D) = \ 



This is a critical exponent. Due to the rearrangement of the Fredholm series, this expansion, as 
similar perturbation expansions, is asymptotic only. 

The relation of parameters A and z to M"^ and e is similar to relations ( p7|) and (|2S 



and 



{D-2) z /•! R{y) Ay 



(40) 



2'kJ-i v^T^V^a {l + zi) 
{D-2) z /■! R{y) dy 



2\iA (1 + z e) 

) = e. (41) 



VI + 2 y (VI + ^ e + Vi + -2 y) 

We will investigate now the positivity of the density function. It is fairly easy to show that in 
the immediate neighborhood of x = —1, where the kernel and inhomogeneity become singular at 
z = 1, the iterated density function is positive definite. Therefore we will concentrate on the region 
of a; + 1 >> I - z. 

The first nontrivial approximation to the Fredholm series is identical to the once iterated solution 
of the integral equation. It provides an asymptotic relation for the range of e. To see this we 
investigate the dependence of A on z. This must be done before one is able to calculate the leading 
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contributions to the density function. If we study solutions at vanishing or smah renormahzed mass. 
M, then the location of the critical point is independent of M in leading orders of the expansion. 

We calculated nonleading order corrections, as well, just to see how much they influence asymp- 
totic relationships. Using ( |40| ) we obtain 

, , / 2 ^ l-(g-2)log[2/(l-z)]/192^2 

^ = ^°Hr^ J 8 + 0.82247 (D-2)/vr^ ^''^ 

When the leading asymptotic part of this relation is substituted into the expression of the density 
function we obtain 



The density function becomes indefinite at large enough D — 2. This is a nonsingular structural 
phase transition, except in the vicinity oi D = 2. Near the transition (D — 2) (log[2/(l — 2)])^ = 



0(1), so at large log[2/(l — z)] the corrections in ( |4^ ) are indeed small. 

Having determined the critical value of z (the maximal z with admissible spectral function) we 
can determine the critical value of e using (|4l|). Substituting into this relation we get 



'log 



^ + - (44) 
3 vr 



Combining (^) and ( ^3|) we get the final expression for the critical value of e 

- 1 / 1024 V2 \ _ 2.331 

~ {D - 2)1/4 27 7r2 ) ~ (Ll- 2)1/4- ^ ^ 

As expected, at = 2 the critical value of e diverges. The prediction of ( pSj ) for the value of ec(-D) 
is not reliable ai D = 4. A better value could be obtained from a thorough numerical study. 

It is interesting to speculate what happens at e > ec{D). It is easy to see that density functions, 
satisfying p{v) > and having finite action still exist. The positive density function that maximizes 
5 at e > Ec may not belong to the class of density functions we have been investigating above. One 
is able to get a hint about the nature of the density function in this region if one considers that 
when e — > edD) — then R(x) develops a zero at some intermediate value of where 
— 1 < Xc < 1. Then at e > edD) the spectrum of eigenvalues should split into two disjunct intervals, 
the density function vanishing at the end points of both intervals. This would be a density function 
that does not belong to the class of functions considered above. Clearly, the change from one to 
the other type of density functions is, an almost probably continuous, structural phase transition. 

To study the nature of these phases it would be very important to derive the dependence of 
large Wilson loops, 

W[X,Y] = (Tr [{a^ f {Uy f (C/.)"^ (Uy)-'']). (46) 

on the area of the loop, XY . This, and the investigation of the model in the strong coupling regime 
of the single link integral are left to a future publication. 
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5 Quenched Theory 



In the above described theory, in the weak couphng phase, the distribution of eigenvalues of the 
gauge matrices is incorrect. To ehminate this problem one must constrain the eigenvalues of the 
gauge field C/^ in a gauge invariant manner. One possible approach is that of Gross and Kitazawa |^ 
that introduces the following modification of the measure: 



(47) 



While this addition to the measure fixes the eigenvalue distribution of matrices C/^, the substitution 
does not affect agreement with the perturbation expansion of the original theory. 
Substituting into the Lagrangian we obtain 



f u 



(48) 



At this point one can introduce the bilinear combinations just as in (^) to get the quenched 
theory for the positive adjoint matrix The action, in terms of the eigenvalues of <I>, (pa takes the 
form 



^ + Nf log (pa - m?(pa - -7r4>l 



+ ^logI(cI>,C^), 



(49) 



where /(<&, Cn) is the Itzykson-Zuber integral [O 



Det 



^2K,(pa COS fcjj 



and where 



A ((^„) A (cos k^) 
Aiaa) = n 



(50) 
(51) 



and is a diagonal matrix formed from the cosines of the momentum variables. 

The density function for the distribution of continuous eigenvalues can proceed the same way 
as it was done in the previous sections. The singular integral equation derived from the variation 
of the action is complicated but can be studied in the weak and strong coupling limits of the 
Itzykson-Zuber integral. 

Finally, we present the expression of Wilson loops in terms of the variables used in this section: 
The Wilson loop of edges X and Y along the x and y axes is defined in the quenched theory as 



{W{X,Y)) 



Tr 
Tr 



(52) 



6 Conclusion 

The investigation of gauge theories induced by scalars in the fundamental representation may 
solve some of the problems that have arisen in theories induced by adjoint representation scalars. 
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Prominent among these is the vanishing of the expectation value of the Wilson loop. This is the 
main reasons why we started an investigation of these theories. 

The most serious difficulty in solving theories induced by fundamental representation scalars has 
been the apparent lack of candidates for a master field. We solved this problems by showing that 
bilinears formed from fundamental representation fields form self adjoint matrices the eigenvalues 
of which are the perfect candidates for a master field. 

Using an Eguchi-Kawai type reduced representation we obtained an integral equation for the 
density of eigenvalues that could be solved exactly in two dimensions at all values of the physical 
parameters. At > 2 we obtained a Fredholm integral equation for the density function. Using 
the perturbative solution of this equation we established that the spectrum undergoes a phase 
transition at a critical value of e = Nf/Nc- Further investigations are needed to examine the 
nature of the phase transition at D > 2 and the size dependence of the expectation value of Wilson 
loops. 

The translation invariant approximation used in most of this paper should be complemented 
by a non-translation invariant approach. Investigation of theories induced by fundamental scalars 
using the Gross-Kitazawa quenching prescription Q is the next step of investigations. 
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